Phase diagram of the anisotropic multichannel Kondo Hamiltonian revisited 
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The phase diagram of the multichannel Kondo Hamiltonian with an XXZ spin-exchange 
anisotropy is revisited, revealing a far richer fixed-point structure than previously appreciated. 
For a spin-^ impurity and k > 2 conduction-electron channels, a second ferromagnetic-like domain 
is found deep in the antiferromagnetic regime. The new domain extends above a (typically large) 
critical longitudinal coupling J* > 0, and is separated from the antiferromagnetic domain by a sec- 
ond Kosterliz-Thouless line. A similar line of stable ferromagnetic-like fixed points with a residual 
isospin- 1 local moment is shown to exist for large J z S> | J± \ > and arbitrary k and s obeying 
\k — 2s | > 1. Here J z is the longitudinal spin-exchange coupling, J± is the transverse coupling, and 
s is the impurity spin. Near the free-impurity fixed-point, spin-exchange anisotropy is a relevant 
perturbation for s > 1/2 and arbitrary k. Depending on the sign of jf — j\ and the parity of 
2s, the system flows either to a conventional Fermi liquid with no residual degeneracy, or to a k- 
channel, spin-^ Kondo effect, or to a line of ferromagnetic-like fixed points with a residual isospin-i 
local moment. These results are obtained through a combination of perturbative renormalization- 
group techniques, Abelian bosonization, a strong-coupling expansion in 1/J Z , and explicit numerical 
renormalization-group calculations. 

PACS numbers: 75.20.Hr, 72.15.Qm 
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INTRODUCTION AND OVERVIEW OF 
RESULTS 



For over the last forty years, the Kondo effect has occu- 
pied a central place in condensed matter physics. While 
earlier studies of the effect have focused on its single- 
channel version realized in dilute magnetic alloys and 
valence-fluctuating systems, later attention has largely 
shifted to its more exotic multichannel variants where de- 
viations from conventional Fermi-liquid behavior can be 
found. The overscreened Kondo effect is a paradigmatic 
example for quantum criticality in quantum-impurity 
systems. Besides the possible relevance to certain heavy 
fermion alloys*^ ballistic metal point contacts)^ scat- 
tering off two-level tunneling systems, 5,6,7 and the charg- 
ing of small metallic grains*^ the overscreened Kondo ef- 
fect is one of the rare examples of interaction-driven non- 
Fermi-liquid behavior that is well understood and well 
characterized theoreticallyiiii The underscreened Kondo 
effect, which might be realized in ultrasmall quantum 
dots with an even number of electrons ^* is a prototype 
for yet another form of unconventional behavior — that 
of a singular Fermi liquid^ 

These vastly different ground states, as well as that 
of an ordinary Fermi liquid, can all be described within 
the single framework of the multichannel Kondo Hamil- 
tonian, which is among the simplest yet richest mod- 
els for strong electronic correlations in condensed-matter 
physics. The multichannel Kondo Hamiltonian of Eq. © 
describes the spin-exchange interaction of a spin-s local 
moment with k identical, independent bands of spin- 
\ electrons. For antiferromagnetic exchange, the low- 
energy physics features a subtle interplay between k and 
s, which could be summarized as follows^ For k = 2s, 



the impurity spin is exactly screened. A spin singlet pro- 
gressively forms below a characteristic Kondo tempera- 
ture Tk, leading to the formation of a local Fermi liq- 
uid. For k > 2s, the impurity spin is overscreened by 
the k conduction-electron channels. The system flows to 
an intermediate-coupling, non-Fermi-liquid fixed point, 
characterized by anomalous thermodynamic and dynam- 
ical properties. The elementary excitations are collective 
in nature, in contrast to the concept of a Fermi liquid. 
For k < 2s, the impurity spin is underscreened. The low- 
energy physics comprises of quasiparticle excitations plus 
a residual moment of size s' = s — k/2. However, it differs 
from a conventional Fermi liquid in the singular energy 
dependence of the scattering phase shift and divergence 
of the quasiparticle density of statesiiSiiiLii Such behav- 
ior was recently termed a singular Fermi liquid. 12 Similar 
qualitative behavior is found for ferromagnetic exchange 
with arbitrary k and s, except that the residual moment 
has the full spin s. 

Inherent to some of the leading scenarios for the real- 
ization of the multichannel Kondo model 5,6,8 is a large 
spin-exchange anisotropy. An XXZ anisotropy is well 
known to be irrelevant both in the single-channel (k = 1, 
s = 1/2) and the two-channel (k = 2, s = 1/2) cases. 
The accepted phase diagram for these two models con- 
sists of an antiferromagnetic and a ferromagnetic do- 
main, separated by a Kosterliz-Thouless line that traces 
J z = — |Jj_| at weak coupling. Here J z and J± are the 
longitudinal and transverse exchange couplings, respec- 
tively. As long as one lies within the confines of the anti- 
ferromagnetic domain, the system flows to the isotropic 
spin-exchange fixed point regardless of how large the ex- 
change anisotropy is. 

Far less explored is the role of spin-exchange anisotropy 
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FIG. 1: (Color online) Phase diagram of the multichannel Kondo model for the different categories of s and k. (a) s = 1/2 
and k > 2. In addition to the well-established ferromagnetic domain for J z < — \J±\, a second ferromagnetic- like domain 
extends to the right and below the red line stretching from (J Z ,J±) = (J*,0). Within bosonization, J* is given by Eq. J37J. 
(b) 1/2 < s < (k — l)/2. In addition to the conventional ferromagnetic domain for J z < — \J±\, the system flows to a line of 
stable ferromagnetic-like fixed points with a residual isospin-i local moment in two opposite limits: (i) for sufficiently small 
Jz > (h) for sufficiently large J z 3> |J±|,D (D being the conduction-electron bandwidth). For weak \J±\ > \J Z \, the 

system flows either to a Fermi- liquid fixed point with a frozen impurity (for integer s), or to an overscreened fe-channel, spin-i 
Kondo effect (for half-integer s). The above weak-coupling behavior is generic to s > 1/2. It extends to all cases described 
below, irrespective of k. (c) s — (k — l)/2. In contrast to (b), the isotropic overscreened fixed point is stable against a weak 
spin-exchange anisotropy (i.e., in the vicinity of the overscreened fixed point), as well as against a sufficiently large J z 3> \ J±\, D. 
(d) s — k/2. The exactly screened fixed point, corresponding to strong coupling, is stable against a sufficiently large anisotropy, 
J z S> | <7j_ |, Z?. However, it is unstable at weak coupling, as described above, (e) s > k/2. All underscreened cases flow to a 
line of stable ferromagnetic-like fixed points with a residual isospin-| local moment both for sufficiently large J z 2> \J±\ and 
for sufficiently small J z > \J±\. 



for either s > 1/2 or i > 2. The stability of 
the overscreened fixed point against weak spin-exchange 
anisotropy has been analyzed in Ref. using confor- 
mal field theory. For k > 1 and either s — 1/2 or s = 
(k — l)/2, the non-Fermi-liquid fixed point was found to 
be stable against a weak spin-exchange anisotropy. This 
has led to the perception that spin-exchange anisotropy 
is irrelevant for these values of k and s. In contrast, spin- 
exchange anisotropy is a relevant perturbation at the 
overscreened fixed point for all other values of 1/2 < s < 
(k — l)/2 (assuming k > 4; for k = 4 it is a marginal per- 
turbation) jiii though the nature of the anisotropic fixed 
points was never explored. Likewise unexplored is the 
effect of spin-exchange anisotropy on the underscreened 
fixed point for s > k/2. 

In this paper, we revisit the phase diagram of the 
multichannel Kondo model with an XXZ spin-exchange 
anisotropy. We find a far more complex picture than 
previously appreciated, including new coupling regimes 
where an XXZ anisotropy substantially alters the low- 
energy physics. Our main findings are as follows. 

• For s — 1/2 and k > 2, a second ferromagnetic- 
like domain is found deep in the antifcrromagnetic 
regime. The new domain extends above a (typically 
large) critical longitudinal coupling J* > 0, and is 
separated from the conventional antiferromagnetic 



(non-Fermi-liquid) domain by a second Kosterliz- 
Thouless line. 

• For spin s > 1/2 and arbitrary k, spin-exchange 
anisotropy is relevant near the free-impurity fixed 
point. For sufficiently small < \pJ±\ < \pJ z \ "C 
1/ks (p is the conduction-electron density of states 
per lattice site), the system flows to a line of sta- 
ble ferromagnetic-like fixed points with a residual 
isospin-i moment. For sufficiently small \pJ z \ < 
\pJ±\ *C 1/ks, the flow is either to a conventional 
Fermi liquid with no residual degeneracy (for half- 
integer s), or to a fc-channel Kondo effect with an 
effective spin-^ local moment (for integer s). Here 
by sufficiently small p J z and pJ± we mean the limit 
J Z ,J± — > with any fixed ratio r — J z /J± ^ ±1- 
The closer \r\ is to one, the smaller the couplings 
must be for these results to apply. 

• For \k — 2s | > 1 and a sufficiently large J z ^> | Jj_| > 
0, the system flows to a line of stable ferromagnetic- 
like fixed points with a residual isospin-i moment. 
Only for k — 2s and k = 2s + 1 are the ex- 
actly screened and overscreened fixed points stable 
against such a large anisotropy, which is marginal 
for k = 2s- 1. 



A compilation of these results is presented in Fig. ^ hi 



the form of phase diagrams for the different categories of 
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s and fc. The phase diagrams for s > 1/2 are still incom- 
plete. There remain extended coupling regimes where 
the low-energy physics is yet to be determined. 

To obtain these results, we employ a combination 
of perturbative renormalization-group (RG) techniques, 
Abelian bosonization, a strong-coupling expansion in 
1/ J z , and explicit numerical renormalization-groupi£ cal- 
culations. Using perturbative RG we first analyze in 
Sec. [n] the limit of weak coupling. In addition to the 
standard RG equations for the dimensionless couplings 
pj z and pJ±, a new Hamiltonian term proportional to 
S z is generated for s > 1/2 and J z =^ J±. Here S z is 
the z component of the impurity spin. Depending on the 
sign of J z — J\, the new Hamiltonian term favors either 
the maximally polarized impurity states (S z — ±s) or 
the minimally polarized ones. This leads to a qualitative 
distinction between J z > \J±\ and |JjJ > \J Z \, and to 
the different low-temperature behaviors described above. 

Proceeding with Abelian bosonization, we next show in 
Sec. IIIII that the anisotropic multichannel Kondo Hamil- 
tonian with s = 1/2 and fc > 1 possesses an exact map- 
ping between the two sets of coupling constants (J z , J±) 
and (J' z , Jj_), where 
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(1) 



The above mapping is restricted to values of J z where 
the right-hand side of Eq. does not exceed ir/2. For 
fc = 2, this condition limits the validity of the mapping 
to J z > 0, in which case Eq. (|TJ simplifies to 



' pJ z ' 



(2) 



For fc > 2, the mapping extends to negative values of J z , 
relating J z > J* = (i/wp) tan(7r/fc) to J m ; n < J' z < and 
vice versa [see Eq. for definition of J m ; n ]. Thus, in 
contrast to common perception, the multichannel Kondo 
Hamiltonian with s = 1/2 and fc > 2 possesses a line of 
stable ferromagnetic-like fixed points for J z > J*. Note, 
however, that J* exceeds the bandwidth for intermediate 
values of fc, and is pushed to weak coupling for fc ^> 1. 

Interestingly, the mapping of Eq. Q is ingrained in 
the Anderson- Yuval approach to the multichannel Kondo 
problem, devised in Refs. andlisl In fact, it was al- 
ready recognized by Fabrizio et al. for fc = 2ji£ but was 
never appreciated to our knowledge for fc > 2. Here we 
provide an explicit operator mapping between the two 
sets of model parameters, for arbitrary fc > 1. 

Since the critical coupling J* predicted by Eq. is 
typically larger than the bandwidth, one might question 
the applicability of bosonization and its unbound linear 
dispersion. Could it be that J* — > oo when the conduc- 
tion electrons are placed on a lattice? To eliminate this 
concern, a strong-coupling expansion in 1/ J z is carried 
out in Sees. II VI and IVl first for s = 1/2 and then for arbi- 
trary s. The strong-coupling expansion not only confirms 
the existence of the new line of stable ferromagnetic-like 



fixed points for s = 1/2 and fc > 2, but further extends 
this result to arbitrary s and fc obeying |fc — 2s| > 1. For 
fc = 2s and fc = 2s + 1, it establishes the irrelevance of 
such a large anisotropy. 

As an explicit demonstration of these ideas, the phase 
diagram of the s = 1/2, fc = 3 model is studied in Sec. I VII 
using Wilson's numerical renormalization-group (NRG) 
method^ The NRG results confirm the phase diagram 
inferred from Eq. (pT|> . including the order of magnitude 
of the critical coupling J* . An extension of the mapping 
of Eq. to a spin-one impurity is presented in turn in 
Appendix lAl We conclude in Sec. IVIll with a discussion 
and summary of our results. 



II. WEAK COUPLING 

We begin our discussion with the limit of weak cou- 
pling, which is treated using perturbative RG. In the 
standard notation, the multichannel Kondo Hamiltonian 
reads 



« = y y y^- h 



Ayr f c t 

2N ^ ^ L fcr 
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s z 
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C-* Ct", 
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(3) 



fc,fc' 



Here, ct creates an electron with wave number fc and 

kna 

spin projection a in the nth conduction-electron channel; 
S is a spin-s operator; J z and J± are the longitudinal 
and transverse Kondo couplings, respectively; and N is 
the number of lattice sites. 

As a first step toward devising a perturbative RG treat- 
ment of the Hamiltonian of Eq. © , we convert the model 
to dimensionless form. To this end, we introduce the 
fermion operators 



1 



^Dp(E)N ^ 



(4) 



which represent the 2fc conduction-electron modes that 
couple to the impurity within the energy shell E = eD. 
Here D is the conduction-electron bandwidth and p{E) is 
the conduction-electron density of states per lattice site: 



(•5) 



The dimensionless operators a\ nG have been normalized 



to obey canonical anticommutation relations: 
ia zna ,a\, n , a ,\ = 5 nn >6 aa i5{E- e'). 



(6) 
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Assuming a box density of states p(E) — pO(D — \E\) 
and omitting all modes that decouple from the impurity, 
Eq. J3J is recast in the form 7i = DTi, where H is the 
dimensionless Hamiltonian 



n=l <7=T,J. 



T 

n=l 
fc 



s% n(7 a ena d£ 



J± \ - 
~9~ 2^ 



cfe / rfe' 



H.c. 



(7) 



Here J 2 = pJ z and Jj_ = p.J± are the dimensionless 
Kondo couplings. 

Focusing on \ J Z \, \ J±\ -C 1/fcs, we treat the Hamilto- 
nian of Eq. JJJ using perturbative RG. The RG trans- 
formation consists of the following three steps. Suppose 
that the bandwidth has already been lowered from its 
initial value D to some value D' = De~ l {I > 0). Further 
lowering the bandwidth to D" — D' (1 — 61) requires the 
elimination of all a Ena degrees of freedom in the interval 
1 — 61 < \e\ < 1. This goal is accomplished using An- 
derson's poor-man's scaling^ At the conclusion of this 
step, all integrations in Eq. Q have been reduced to the 
range —(1 — 61) < x < 1 — 61. The RG transformation 
is completed by (i) rescaling H — > H/(l — 61) to account 
for the reduced bandwidth, and (ii) restoring the original 
integration ranges in Eq. JIJ) by converting to 

a e n<j — > a ena = (1 — £J) a (i-<5i)eno- (8) 
This latter step allows us to write 

/l-Sl pi 
aenade = (1 - 61) 1/2 / a ena de. (9) 
-(1-51) J-l 

The dimensionless Hamiltonian is recast in this manner 
in a self-similar form, but with renormalized couplings. 
These obey a set of coupled differential equations, speci- 
fied below. 

For s = 1/2 we recover the familiar RG equations 



dl 
dJ ± 
~df 



T 2 



J±Jz 



(10) 
(11) 



showing that spin-exchange anisotropy is irrelevant for 
weak J z > — \J±\, irrespective of the number of chan- 
nels k. A different qualitative picture emerges for spin s 
larger than one-half. Here a new Hamiltonian term AS 2 
is generated within TL. Starting from zero, A renormal- 
izes according to 



which supplements Eqs. IjlOfl and (|11|) for J z and J±. 

Restricting attention to s > 1/2, we now analyze the 
ramifications of the new coupling A as a function of s, fc, 
and the bare Kondo couplings J z and Jj_. 

Since J 2 — J\ is conserved under the RG, a negative 
(positive) A is generated if the bare Kondo couplings 
satisfy \J Z \ > \J±\ (\J Z \ < \J±\). For a given ratio 
Jz/J± 7^ ±1 and sufficiently small |Jj_|, the coupling 
|A| approaches unity well before any Kondo temperature 
can be reached. This has the effect of freezing all but the 
lowest-lying spin states. For \J Z \ > |Jj_| (negative A), 
only the maximally polarized states S z = ±s are thus 
left. For | J z \ < \J±\ (positive A), the picture depends on 
s: for half-integer s, the two degenerate states S z = ±1/2 
are selected; for integer s, only the state S z = remains. 

Depending on which case is realized, a different effec- 
tive low-energy Hamiltonian emerges. Consider first the 
case \J Z \ > \J±\ (negative A). Introducing the isospin 
operators 



l*X*l 



*)(-*!, 



±s)( T s\ 



(13) 



(14) 



the resulting low-energy Hamiltonian contains the term 



H, 



s j z y. 



de I ds' 



which follows from projection of the spin-exchange in- 
teraction of Eq. JJJ) onto the S z — ±s subspace. Here 
J z is the running coupling constant at the scale where 
|A| approaches one. In contrast to Eq. l(T5|) . terms that 
flip the isospin f involve the creation and annihilation 
of at least 2s conduction electrons. This follows from 
the fact that a flip in f changes S z by ±2s. Since the 
Hamiltonian conserves the total spin projection of the 
entire system (impurity plus conduction electrons), such 
a process must be accompanied by an opposite spin flip 
of 2s conduction electrons in the vicinity of the impurity. 
As result, isospin-flip terms have the scaling dimension 
Aj^ > 2s > 2 about the free-impurity fixed point, ren- 
dering them irrelevant. The resulting fixed-point struc- 
ture corresponds then to a line of ferromagnetic-like fixed 
points with a residual isospin- i local moment, character- 
ized by a finite J z . Note that this picture is insensitive 
to the sign of the Kondo couplings. It equally applies to 
positive and negative J z and J±. 

Proceeding with the case \J Z \ < \J±\ (positive A), we 
first consider half-integer s. Here the two states selected 
by A are S z = ±1/2. Similar to Eqs. and ifHjl. we 
introduce the isospin operators 

r l|l/2)<l/2|-i|-l/2)<-l/2|, (16) 



dA 
~dT 



A - feln(2) 



t2 _ t2 



(12) 



±1/2)(±T/2|, 



(17) 
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which now relate to the states \S Z = ±1/2). Projection of 
Eq. J7J) onto the S z = ±1/2 subspace yields the following 
effective spin-exchange interaction: 



n, 



J, 



E 



de / de' 



J± 

2 



7^E 



i ,i 
de / de' 
-l J-i 



H.c. 



(18) 



with 7 = y/s(s + 1) + 1/4 > 1. Once again, J z and Jj_ in 
Eq. (|18fl are the running coupling constants at the scale 
where A approaches one. 

Equation (|18fl has the form of a fc-channel Kondo 
Hamiltonian with an effective isospin-^ local moment. 

Since | J z \ < -f\J± \ <C 1, one lies in the confines of the an- 
tiferromagnetic domain. Hence, irrespective of the orig- 
inal spin s, the system flows to the overscreened fixed 
point of the fc-channel Kondo effect with spin s' = 1/2. 
(For fc = 1, the flow is to the strong-coupling fixed point 
of the conventional one-channel Kondo effect). Exclud- 
ing the case s = fc/2 — 1/2 with antifcrromagnetic ex- 
change, the resulting low-energy physics differs markedly 
from that of the isotropic model, whether ferromagnetic 
or antiferromagnetic. Most strikingly, the underscreened 
fixed point for J z = J± > and s > fc/2 > 1/2 gives 
way to an overscreened fixed point for any given ratio 
— 1 < J z /Jj_ < 1 and sufficiently small |JjJ. 

Of the different possible cases for A and s, the simplest 
picture is recovered for \J Z \ < \J±\ and integer s. Here 
the impurity spin is frozen in the S z = state, loosing 
its dynamics. This results in a conventional Fermi-liquid 
fixed point with neither non-Fcrmi-liquid characteristics 
nor a residual local-moment degeneracy. 



III. EXACT MAPPING FOR s = 1/2 

As is evident from the discussion above, there is a qual- 
itative difference between s = 1/2 and s > 1/2 con- 
cerning the effect of spin-exchange anisotropy for weak 
coupling. In this section we focus on s — 1/2 and use 
Abelian bosonization to derive the mapping of Eq. 
Our starting point is the continuum-limit representation 
of the multichannel Kondo Hamiltonian 



fe />oo 

H = E i ^ iVF / ^} ltT ( x )9x^na(x)dx 
1 _ * , J— oo 



(19) 



71 = 1<T = T,X 

fc 



n=l 
k 



+ Et" ^i(o)Vm(o)s+ + S-< t (0)iM0) 



n=l 

9 eh 



2 

k 

E 

71=1 



Vd T (o)VM(o)-v4i(o)VM(o) 



9ihS z 



dx, 



written in terms of the left-moving one-dimensional fields 
ipta( x ) with n — 1, • ■ • , fc and a =|, J.. Here, S is a spin-i 
operator; a is a short-distance cutoff corresponding to a 
lattice spacing; a; is a fictitious coordinate conjugate to 
fc = £7r/a; and h is an applied magnetic field. For the 
sake of generality, we allow for different impurity and 
conduction-electron Lande g-factors, <?i and g e , respec- 
tively. 

To treat the Hamiltonian of Eq. i|19[l . we resort to 
Abelian bosonization. According to standard prescrip- 
tions^ 2fc boson fields are introduced — one boson field 
*&na{x) for each left-moving fermion field ip na {x). The 
fermion fields are written as 



1 



V2t 



(20) 



where the fields obey 

[®ncr(x),$n'a'(y)] = -i6 nn > 6^' it sign(x - y). (21) 

The ultraviolet momentum cutoff a -1 = n/a is related 
to the conduction-electron bandwidth D and the density 
of states per lattice site p through D = hvp/a and p — 
1/(2D) — a/(2hvp), respectively. The operators P na 
in Eq. I|20|l are phase-factor operators, which come to 
ensure that the different fermion species anticommutc. 
Our explicit choices for these operators are 



(22) 



where Nj a is the number operator for electrons in channel 
j with spin projection a . 

In terms of the boson fields, the multichannel Kondo 
Hamiltonian takes the form 



n=l cr=T,i 
fc 



{\7§ na {x)Ydx 
^ ^ {e 4 [*"T(0)-*M(0)] 5 - +H .c.} 



k 

s *-r- E I v$ ™t(o) - v$ ni (o)] S z - 9l hS z 



ir*p 



9 eh X - 
4tt <H 



[V$„ T (a;) - V$ ni (cc)]dx, (23) 



where 



arctan 



irpj z 



(24) 



is the parallel-spin phase shift induced by J z in the 
absence of J±. Note that S z is bounded in magni- 
tude by 7r/2, which stems from the cutoff scheme used 
in bosonization. Although the bosonic Hamiltonian of 
Eq. does support larger values of \S Z \, this param- 
eter must not exceed ir/2 in order for 7i to possess a 



6 



fermionic counterpart of the form of Eq. I|19[) . We shall 
return to this important point later on. 

At this stage we manipulate the bosonic Hamiltonian 
of Eq. i|23|) through a sequence of steps. We begin by 
converting to 2k new boson fields, consisting of 



$ s (x) = -L= [*n T (*) - *M (x)] (25) 



plus 2k — 1 orthogonal fields: with fi= 1, ■ ■ ■ ,2k — 

1. The orthogonal fields $/x(a;) are formally expressed as 



n=l <t=T,1 

where the real coefficients e^ CT obey 

A; 

n=l cr=Ta 



(26) 



e nT nj. 



= 0. 



(27) 



(28) 



The precise form of the coefficients e^ CT is of no practical 
importance to our discussion, and need not concern us. 
Their choice is not unique. In terms of the new fields, 
the combinations $ n t(a;) — $ n i (a;) take the form 



-$ s (x)+ip n (x), (29) 



where (^ rl (a;) is some linear combination of the fields 
&fj,(x) with /i = 1, • • • , 2k — 1. Once again, the precise 
form of the <p n (x) combinations is of no real significance 
to our discussion. We shall only rely on them being or- 
thogonal to $ s (ab- 
using the new boson fields defined above, the Hamil- 
tonian of Eq. (|23[l is converted to 



H 



Tivf 
4tt 



2k-l 



(v$ s ) 2 + ]T (v*„ 



dx 



+ S 2 



, . 2W$ S (0)S, - gihS, 



4tt 



2k / V<P s (x)dx. 



(30) 



Next thecanonical transformation Ti' = UHTP with U = 
exp 



H/ I^s(0)S z is applied to obtain 



w = 



Tivf 

Air 



2fc-l 



(/.?: 



+ £ J -L { e -^a(0)+^»(0) 5 - + Rc .| 



^ 2 

71=1 



^ -^-V2fcV$ s (0)5, - ( 9i - 2g e )hS z 



9 —^2~k / V$,(s)da;. 



(31) 



Here we have omitted a constant term from H', and made 
use of the identity a/irp — 2ft, in writing the term 
V$ s (0)S z . Finally, the transformation $> s (x) — > —& a (x) 
is introduced, which yields 



H' 



Tivp 

47T 

fc 



2k-l 

■E 



£^{ e VfMo)+Wo)< r+ILc .} 

n— 1 

/ 2fcV$ s (0)5 z + (2 9e -^)^ 



47T 



2fc / V$ s (x)dx 



(32) 



with <^ = n/k — 5 Z . 

The Hamiltonian of Eq. (|32|l is identical to that of 
Eq. (|30(l , apart from a renormalization of certain param- 
eters: 5 Z — ► (5^,, ft, — * —ft, and — * 2(7 e — <?j. As long 
as |^| < 7r/2, one can revert the series of steps lead- 
ing to Eq. 13Ufl . to recast the Hamiltonian of Eq. IMl'l) 
in fermionic form. The end result is just the original 
multichannel Kondo Hamiltonian of Eq. (|f 9J) with the 
following renormalized parameters: 



7^7' 



■ tan 



7T/9 



f- 

Vfc 



ft -> 2g e - gi 



(33) 



(34) 



(35) 



For zero magnetic field, this establishes the mapping 
of Eq. (JTJ, including the restriction to values of J z where 
the right-hand side of Eq. 0) does not exceed ir/2. The 
latter condition is just a restatement of the requirement 
\8'z\ 7r /2- We now analyze in detail the ramifications of 
this restriction as a function of the number of channels 
k. 

For k — 1 (single-channel case), S' z exceeds ir/2 for all 
— 00 < J z < 00. Thus, the mapping of Eq. (JTJ does not 
apply to the single-channel Kondo effect, in accordance 
with known results. For k — 2 (two-channel case), the 
required condition is met for all J z > 0, mapping weak to 
strong coupling and vice versa [see Eq. (J5J]. Hence, the 
mapping of Eq. can be viewed as an anisotropic vari- 
ant of the weak-to-strong-coupling duality of Nozieres 
and BlandiniSi 



7 



The most interesting case occurs for k > 2, when 
Eq. 0} extends to all J z > J min with 



(36) 



Jmin = tan - - - < 0. 

np V 2 k 1 

In particular, the range J z > J* with 

4 



JZ 



■ tan 



Tip 



f- 



> 



(37) 



is mapped onto the negative-coupling regime J m i„ < 
J' z < and vice verse. Thus, the Kosterliz-Thouless line 
separating the antiferromagnetic and ferromagnetic do- 
mains is duplicated from J z = — | JjJ to J z = J* +Ck\ J±\ 
with 



Cfc 



1 



1 



l + tan 2 (7r/fc) l + ^pJ*/^) 2 



(38) 



Here we have assumed p\J z — J*\,p\J±\ <C 1, in order 
for the weak-coupling parameterization of the Kosterliz- 
Thouless line to apply. 

The resulting phase diagram for k > 2 is plot- 
ted in Fig. 12a). For J z > J* + C k \J±\, the mul- 
tichannel Kondo Hamiltonian flows to a line of stable 
ferromagnetic-like fixed points, rendering the non-Fermi- 
liquid fixed point of the model unstable against a large 
enough anisotropy. This behavior contradicts the com- 
mon perception of spin-exchange anisotropy as irrelevant 
for s = 1/2. Note that the same phase diagram emerges 
from the renormalization-group equations derived by Ye 
using the equivalent Anderson- Yuval approach. 18 How- 
ever, the newly found line of stable ferromagnetic-like 
fixed points has eluded Ye. 



IV. STRONG-COUPLING EXPANSION FOR 

s = 1/2 

A potential concern with the above picture for s = 1/2 
has to do with the validity of the bosonization approach 
used. Since J* exceeds the bandwidth D for intermediate 
values of k, one might wonder to what extent is bosoniza- 
tion (or the Anderson- Yuval approach for that matter) 
justified for such strong coupling. In fact, the very usage 
of the continuum-limit Hamiltonian of Eq. I|19|) with its 
unbounded linear dispersion can be called into question. 
Could it be that J* — > oo when the conduction electrons 
are placed on a lattice? 

Not withstanding the observation that J* for k ^> 1 
is pushed to weak coupling, to firmly establish the phase 
diagram of Fig. ^a) one must show that the new line of 
stable ferromagnetic-like fixed points extends to proper 
lattice models for the underlying conduction bands. This 
is the objective of the following analysis, which focuses on 
the limit of a large longitudinal coupling, J z ^> D, |Jj_|- 

As a generic lattice model for the conduction bands, we 
consider a spin-i impurity moment coupled to the open 



end of a semi-infinite tight-binding chain with k identical 
conduction-electron species: 



n 



UXJ ft, 

EE E \_ e jfjnafjna- 

3=0 r»=l<7=T4 

"K? \f(j+l)naf 

k 

^E^i/onT^ + H.c." 

71= 1 
J k 
71=1 

j—0 71=1 



H.c. 



9ihS z 



(39) 



Here ej and tj, respectively, are the on-site energies and 
hopping matrix elements along the chain. Any lattice 
model with identical noninteracting bands can be cast 
in the form of Eq. I|39|) using a Wilson-type construc- 
tion^ Different lattice models are distinguished by the 
tight-binding parameters tj and tj, the largest of which 
determines the bandwidth D. For example, particle-hole 
symmetry demands that ej — for all sites along the 
chain. In the following we assume a large longitudinal 
coupling, J z 3> D, |J_l|, and expand in powers of 1/J 2 
to derive an effective low-energy Hamiltonian at energies 
far below J z . 

For ksT <C J z , the fermionic degrees of freedom at site 
zero bind tightly to the impurity so as to minimize the 
J z interaction term. There are two degenerate ground 
states of this interaction term: 



ii/u=t) 



(40) 

corresponding to the total spin projections 6J otaJ = ±(k— 
l)/2. All excited eigenstates of the J z interaction term 
are thermally inaccessible, being removed in energy by 
integer multiples of J z /4. Defining a new isospin operator 
f according to 



(41) 



p=± 



[compare with Eqs. ©-lO and Eqs. (O-JUJ], the ef- 
fective low-energy Hamiltonian takes the form 7i c ham + 
H 



mag 



Hint, where 7i c hain is the tight-binding Hamil- 
tonian of the truncated chain with site j = removed, 
H mag is the magnetic-field term 



Wmag — 9t 



oo k 



j = l 71=1 



(42) 

with g T = kg e — gi, and 7ii n t contains all finite-order 
corrections in either 1/J 2 or J± (or both). 



8 



The explicit form of the Hamiltonian term Hint de- 
pends on the number of conduction-electron channels k. 
For k = 1 and up to linear order in 1/J Z , it takes the 
form 



Hint = Jj_T x + ~^ T z 



/it /it - fl\ hi 



where A z equals 



4/ 



< D. 



(43) 



(44) 



Here we have omitted the redundant channel index 



from within f\ nrT , i.e., we have set f\ na 



/ lCT . Since the 

A z term is exactly marginal, the low-temperature physics 
is governed for ft, = by the spin- flip term Jj_, which 
favors the singlet state [|+) — |— )] for J± > 0. In this 
manner one recovers the characteristic spin singlet of the 
single-channel Kondo effect, restoring thereby SU(2) spin 
symmetry as T — + 0. A local magnetic held breaks the 
emerging SU(2) spin symmetry as it physically should by 
introducing weak spin-dependent scattering at the open 
end of the truncated chain. 

The situation is somewhat different for k — 2. In this 
case Hint acquires the form 



Ht 



A, 



Tz E [/l«T/ lrl T — /ln|/ 
2 



(45) 



where \ z is still given to order 1/J Z by Eq. I|44|) . and Aj_ 

scales as 



A, 



J* ' 



(46) 



Here we have omitted higher order interaction terms 
within Hint, and restricted attention to particle- hole sym- 
metry. Away from particle-hole symmetry, an additional 
potential-scattering term of the form Aei J2 n finer hna is 
generated at second order in 1/ J z . 

Equation (|45|l has the same exact form as the origi- 
nal spin-exchange interaction in Eq. (|39|) . but with two 
modifications: (i) The coupling constants J z and Jj_ have 
been pushed to weak coupling; (ii) The site j = has 
been replaced with j = 1. The physical content of the 
local moment f has also changed. Importantly, since 
D ^> X z ^> | Aj_ | one lies within the confines of the antifer- 
romagnetic domain, rendering Hint a relevant perturba- 
tion. Hence, in accordance with the bosonization treat- 
ment of Sec. IIHI the strong-coupling limit J z 3> D maps 
onto weak coupling, extending the duality of Nozieres 
and BlandinSi to large spin-exchange anisotropy. More- 
over, since to ~ 1/p for conventional lattice models, then 
A 2 of Eq. ijllj) is comparable to J z of Eq. @ . The main 
difference as compared to bosonization pertains to the 



transverse coupling J±, which renormalizes to Aj^ cx 1/ J z 
in the strong-coupling expansion, but is left unchanged 
within bosonization. Excluding this rather minor dis- 
crepancy, the two approaches are in close agreement with 
one another. 

The crucial difference for k > 2 has to do with the 
dynamic components of Hint , responsible for flipping the 
isospin t. To see this we note that the total spin pro- 
jections S* otal of the states |+) and |— ) differ by k — 1. 
Since the Hamiltonian of Eq. (|39fl preserves the overall 
spin projection of the entire system, then the flipping of 
t z from plus to minus or vice versa must be accompa- 
nied by an opposite spin flip of k — I electrons along the 
truncated chain. Hence the leading dynamic term in Hint 
only shows up at order (l/Jz) 2 ^" -1 ), taking the form 



E II fl]hni + a.c. 



with 



A 



J± 



2(k-l) 



(47) 



(48) 



In contrast to the dynamic part, the leading static com- 
ponent of Hint remains given by the A z term of Eq. I|45|l , 
except for the summation over n which now runs over all 
k channels. To linear order in 1/J Z the coupling A z is in- 
dependent of k, being given by Eq. (|44() . Once again, an 
additional potential-scattering term is generated at order 
(1/J Z ) 2 away from particle-hole symmetry. 

Since the Aj^ term of Eq. H47|) involves the creation 
and annihilation of k — I electrons at the open end of the 
truncated chain, it has the scaling dimension Aj^ — fc— I 
with respect to the J z — » oo "free" Hamiltonian. Hence 
this term is irrelevant for k > 2. The same holds true 
of all higher order dynamical terms generated, as these 
likewise contain at least k — 1 creation and k — I annihila- 
tion operators of electrons localized along the truncated 
chain. Similar to the case s > 1/2 and \J Z \ > \J±\ in 
Sec. EI the resulting fixed-point Hamiltonian for h = 
corresponds then to a finite longitudinal coupling A z 
but zero A^, in perfect agreement with the results of 
bosonization. In fact, an identical scaling dimension 
Aj^ — > A; — 1 is obtained in the Anderson- Yuval approach 
for pj z 3> lr^ reinforcing the qualitative agreement be- 
tween the two approaches. Most importantly, a line of 
stable ferromagnetic-like fixed points is seen to exist for 
any lattice model with sufficiently large J z > 0, just as 
predicted by bosonization. 

Evidently, the strong-coupling expansion confirms the 
results of bosonization for all three cases: k = 1, k = 2, 
and k > 2. Moreover, it provides a transparent physical 
picture for the source of distinction between the three 
cases. It all boils down to the nature of the local spin 
configurations selected by a sufficiently large J z > 0. We 
therefore conclude that Fig. ^a) correctly describes the 
phase diagram of the spin-^ multichannel Kondo model 
with k > 2 conduction-electron channels. 
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V. STRONG-COUPLING EXPANSION FOR 
ARBITRARY SPIN s 

Our discussion in the previous two sections was con- 
fined to a spin-i impurity. An appealing feature of the 
strong-coupling expansion in 1/J Z is that it can easily be 
generalized to arbitrary spin s. This is the goal of the 
present section. 

The basic considerations for arbitrary s > are quite 
similar to those for s = 1/2. As before, the J z interaction 
term possesses two degenerate ground states. We label 
these states according to 



k 



n /onTi^* 



for fc > 2s, and 

k 

i+>=rWmi^ 



(49) 



OnT 1 



(50) 

for fc < 2s. With this convention, |±) has the total spin 
projection S* otal = ±\k/2 — s\. Defining the isospin f 
according to Eq. (|41|l . the effective low-energy Hamilto- 
nian is written as W c ham + H mag + H- m t, where H c hain 
is the tight-binding Hamiltonian of the truncated chain 
with site j = removed, 7i m ag is the magnetic-field term 
of Eq. (|42|l . and 7i; n t contains all finite-order corrections 
in either 1/J Z or (or both). The sole modification 
to W ma g of Eq. is in the effective g-factor g Tl which 
changes from g e k — gi for s = 1/2 to ±(g e k — 2sgi) for 
arbitrary s. Here the plus (minus) sign corresponds to 
k > 2s (fc < 2s). 

Similar to the case s = 1/2, the delicate interplay be- 
tween fc and s enters through the Hamiltonian term Hint ■ 
Let us separate the discussion of the dynamic and static 
components of 7ii n t, as these depend differently on s and 
fc. The leading static component of Hint remains given by 
the A z term of Eq. I|45[) . except for the summation over 
n which now runs over all fc conduction-electron chan- 
nels. The coupling A z does depend on fc — 2s, however 
only through its sign. For fc > 2s it is given to order 
1/J Z by Eq. I|44|l . corresponding to an antiferromagnetic 
interaction. For fc < 2s it acquires an additional mi- 
nus sign, corresponding to a ferromagnetic interaction. 
Apart from its overall sign, the leading static component 
of 'Hint is independent of s. 

Moving on to the dynamic part of H- m t, its leading- 
order component displays a more elaborate dependence 
on fc and s. As noted above, the total spin projections of 
the states |+) and |— ) differ by |fc — 2s\. Consequently, 
the flipping of r 2 from up to down or vice versa must be 
accompanied by an opposite spin flip of |fc — 2s| electrons 
along the truncated chain, else the total spin projection of 
the system is not conserved. This consideration dictates 
the following form for the leading dynamical term: 



where 0^ s = (O ks Y is a channel-symmetric operator 
that creates |fc — 2s | spin-up electrons and annihilates 
|fc — 2s| spin-down electrons at, or close as possible to, 
the open end of the truncated chain. 

The operator Of is generally too complicated to write 

down. It greatly simplifies in two cases. For fc = 2s, O ks 
reduces to the unity operator; For s = 1/2 and fc > 1, 
it is given by Eq. (|47|l . As for the coupling A_l , it scales 
differently for fc > s and fc < s. For fc > s, Aj^ behaves 
as 



J± 



J± 

J, 



2|fc-2s| 



(52) 



Ax 



01 



(51) 



For fc < s, it involves a higher power of 1/J Z , as electrons 
farther into the chain must participate in the flipping of 
t z . With the possible exception of fc = 2s, the transverse 
coupling | Ax | is parametrically smaller than |A Z |, a fact 
that will have important implications later on. 

As a function of s and fc, Ax has the scaling dimension 
A± = \k — 2s | with respect to the J z — » oo "free" Hamil- 
tonian. For h — 0, this yields the following classification 
of the low-energy physics. 

(i) s = k/2 — In this exactly screened case, Ax acts 
as a local transverse magnetic field, lifting the two-fold 
degeneracy of |±). Although the local state favored by 
Ax is generally not an SU(2) spin singlet, the low-energy 
physics is identical to that for isotropic antiferromagnetic 
exchange, as can be seen from the boundary condition im- 
posed on the truncated chain. Specifically, a local Fermi 
liquid progressively forms below the Kondo temperature 
T K ~ Ax- 

(ii) s = k/2 — 1/2 — In this overscreened case, 
the Hamiltonian term Hi n t has the same exact form as 
Eq. 1)45(1 . except for the summation over n which now 
runs over all fc conduction-electron channels. Hence the 
system is described by a weakly coupled fc-channel Kondo 
Hamiltonian with A 2 ^> |Aj_| > 0. The effect of the 
large spin-exchange anisotropy in the original Hamilto- 
nian is to reduce the effective impurity moment from s 
to 1/2. Based on the perturbative RG analysis of Sec.lTTl 
the resulting Hamiltonian flows to the overscreened fixed 
point of the fc-channel, spin-i Kondo Hamiltonian, which 
is equivalent in turn to that of the fc-channel Kondo 
model with s = k/2 — 1/2. As in the exactly screened 
case discussed above, a large spin-exchange anisotropy 
J z ^> D, | Jj_| is seen to be irrelevant for s = k/2 ~ 1/2. 

(iii) s = k/2 + 1/2 — Similar to the previous case, 
the system is described by a weakly coupled fc-channel 
Kondo Hamiltonian with s — > 1/2. However, the effective 
coupling is now ferromagnetic: — A z 3> |Ax| > 0. Conse- 
quently, the flow is to a line of stable ferromagnetic-like 
fixed points with finite A z but zero Ax • The resulting low- 
energy physics is that of a singular Fermi liquid 1 - plus 
a residual isospin r. It differs from that of an isotropic 
spin-exchange interaction only in the residual A z inter- 
action. Since the latter term is marginal, so is the large 
spin-exchange anisotropy for this underscreened case. 
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FIG. 2: (Color online) Phase diagram of the multichannel 
Kondo model with J z 3> \J±\,D, as a function of k and s. 
For 2s = k and 2s = k — 1, the low-energy physics is sta- 
ble against such a large spin-exchange anisotropy. The fixed 
point remains that of an exactly screened (2s = k) or an over- 
screened (2s = k — 1) Kondo effect. For \k — 2s| > 1, the low- 
energy physics is unstable against such a large spin-exchange 
anisotropy. The system flows to a line of stable ferromagnetic- 
like fixed points with a residual isospin-i local moment. For 
2s = k + 1, a large spin-exchange anisotropy is marginal. The 
system still flows to a line of stable ferromagnetic-like fixed 
points, however these differ from the isotropic underscreened 
fixed point only by the marginal operator A z (see text). 



(iv) | k — 2s | > 1 — In this effectively underscreened 
case, the scaling dimension Aj_ exceeds one. Hence Aj^ is 
irrelevant, as are all higher order dynamical terms gener- 
ated within Hint ■ The system thus flows to a line of stable 
ferromagnetic-like fixed points, characterized by a finite 
A 2 but zero Aj_ . The low-energy physics is again that of a 
(potentially singular) Fermi liquid plus a residual isospin 
r. For s < k/2— 1/2, this differs markedly from the over- 
screened fixed point of the isotropic spin-exchange model. 
Also for s > k/2 + 1/2 this differs from the underscreened 
fixed point of the isotropic model, as the residual local- 
moment degeneracy is two instead of 2s — fc + 1 > 2. 
Importantly, in all cases the resulting low-energy physics 
is insensitive to the sign of k — 2s, in stark contrast to 
the isotropic case. A large spin-exchange anisotropy of 
the form J z ^S> D,\J±\ is therefore relevant with respect 
to the isotropic fixed point for all \k — 2s\ > 1. 

The phase diagram of the multichannel Kondo Hamil- 
tonian with J z 3> | J±\,D is summarized in Fig. [21 as a 
function of k and s. Excluding the exactly screened line 
2s = k and the overscreened line 2s = k — 1, the model 
flows to a line of stable ferromagnetic-like fixed points 
with a residual isospin-^ local moment for all k and s. 
For 2s = fc+1, the ferromagnetic-like fixed points and the 
isotropic underscreened fixed point are equivalent. They 
only differ by the marginal operator X z . For 2s > k + 1, 
the ferromagnetic-like fixed points and the isotropic un- 



derscreened fixed point are distinctly different, possessing 
a different residual degeneracy. 



VI. NRG STUDY OF s = 1/2, k = 3 

Although the strong-coupling expansion unequivocally 
confirms the existence of a line of stable ferromagnetic- 
like fixed points for large J z and \k — 2s| > 1, it can- 
not access the entire phase diagram of the anisotropic 
multichannel Kondo Hamiltonian. In particular, the sec- 
ond Kosterliz-Thouless line predicted by bosonization for 
s = 1/2 and k > 2 [see Fig.QJa)] lies beyond the scope of 
this approach. In this section, we use Wilson's numeri- 
cal renormalization-group (NRG) methodic to conduct a 
systematic study the phase diagram of the anisotropic 
multichannel Kondo Hamiltonian with s = 1/2 and 
k = 3. 

In conventional formulations of the NRG^ one consid- 
ers a particular choice for the tight-binding parameters 
in Eq. (JSSJ, given by Cj = and t 3 = DaA^^ 2 ^ with 



Da = ya + A- 1 ), 



1 - a-(j +i ) 

v /(l-A-( 2 J+ 1 ))(l- A-( 2 J+ 3 ))' 



(53) 



(54) 



Here A > 1 is a discretization parameter. For A — * 1 + , 
the resulting model describes an impurity spin locally 
coupled to k identical conduction bands with a symmetric 
box density of states p(E) = (l/2D)8{D - \E\). For 
A > 1, the model represents an impurity spin coupled to 
a logarithmically discretized version of the same bands*^ 
The first step in the NRG approach is to recast the 
Hamiltonian TL as the limit of a sequence of dimensionless 
Hamiltonians H.n'- 



H= lim {DaA-^-^Hn} 



with 
TLn 



A (JV-l)/2 



J± 



2D A 



X] (/Li/onT 



H.c. 



n=l 



(55) 



(56) 



+ 



J, 



k 

^s z 



n— 1 

7V-1 fe 



3=0 n=l <7=T,i 



The finite-size Hamiltonians H.n are then diagonalized 
iteratively using the NRG transformation 



H N+1 = VAHn + E E ^ {fU+DnjNna + H.C.} . 

(57) 



ra=l<r=T,J. 
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Here the prefactor A^ -1 )/ 2 in Eq. i|56[l guarantees that 
the low-energy excitations of Hn are of order one for all 
N. The approach to a fixed-point Hamiltonian is signaled 
by a limit cycle of the NRG transformation, with Hn+2 
and Hn sharing the same low-energy spectrum. 

It practice, it is impossible to keep track of the ex- 
ponential increase in the size of the Hilbert space as a 
function of the chain length N. Hence only the lowest 
N s eigenstates of Hn are retained at the conclusion of 
each iteration^ The retained states are used in turn to 
construct the eigenstates of Hn+i- The truncation er- 
ror involved can be systematically controlled by varying 
N s . However, since the size of the Hilbert space increases 
by a factor of 2 2fe with each additional iteration, only 
a moderately small number of conduction-electron chan- 
nels can be treated reliably using present-day computers, 
typically no more than 3 or 4 channels. In the following 
we focus on k — 3 and s = 1/2, which is the simplest 
version of the multichannel Kondo Hamiltonian that is 
both amendable to the NRG and predicted to display 
the phase diagram of Fig. ^a). 

To cope with the large computational effort involved in 
exploring the phase diagram of the three-channel Kondo 
model with spin-exchange anisotropy, we used rather 
large values of A (either A = 3orA = 5), and applied 
an alternative truncation scheme to the one customarily 
used. Since each of the NRG Hamiltonians Hn is block 
diagonal in the conserved quantum numbers^ the com- 
putational time is mostly governed by the largest block to 
be diagonalized. Thus, instead of retaining a fixed num- 
ber of eigenstates of Hn, at each iteration we adjusted 
the threshold energy for truncation so that the largest 
block did not exceed A^iock states* 2 - 3 For A = 3 we used 
A^biock — 100, while for A = 5 we set iVbi oc k = 80. This 
approach gave rise to variations in the number of states 
retained at the conclusion of each iteration. Typically, 



3000 to 4000 states were retained. In some extreme cases 
down to 2500 states and up to 5500 states were kept. We 
have verified that the threshold energies selected in this 
way were sufficiently high so as not to spoil the accuracy 
of the calculations. 

Figure depicts the finite-size spectra obtained for 
pj± = 0.008 and different values of pj z > 0. Here 
A = 3 was used. Up to a critical coupling pJ z sa 1.26, 
the system always flows to the non-Fermi-liquid fixed 
point of the three-channel Kondo effect. Indeed, the 
fixed-point spectrum is independent of J z < J z , and is 
identical to that for isotropic spin-exchange interaction 
(left-most panel). There is excellent agreement with the 
finite-size spectrum of the three-channel Kondo effect ob- 
tained from conformal field theory, 25 whose energy levels 
are marked by arrows in Fig. [3J Here the slight split- 
ting of NRG levels near the dimensionless energy 1.6 is 
a A-dependent feature. This splitting is reduced upon 
decreasing A, and should disappear for A — > 1 + . 

A different picture is recovered for J z > JJ. (i) As 
demonstrated in Fig. for pj z = 1.36 and pJ z = 6.4, 
the finite-size spectrum remains Fermi-liquid-like down 
to the lowest energies reached (E ~ 1CK 120 -D in some 
extreme runs), without crossing over to the overscreened 
fixed point of the three-channel Kondo effect, (ii) The 
finite-size spectrum varies continuously with J z > J°, 
suggesting the flow to a line of fixed points connected 
by a marginal operator, (iii) For J z 3> J z , the quan- 
tum numbers and degeneracies of the NRG levels are in 
excellent agreement with those anticipated based on the 
strong-coupling analysis of Sec. IIVI (iv) The NRG level 
flow confirms that channel asymmetry is a marginal per- 
turbation for J z > J c z (not shown), in stark contrast to 
the regime J z < J z . These features are all consistent 
with the flow to a line of stable ferromagnetic-like fixed 
points as predicted by bosonization. 




pj = 1.36 



pj = 6.4 



FIG. 3: NRG level flow (odd iterations) for s — 1/2, k = 3, pJ± = 0.008, and increasing values of pJ z . Here A = 3 was 
used. For pj z < pj z ~ 1.26 (left two panels), the system flows to the non- Fermi- liquid fixed point of the three-channel Kondo 
effect, irrespective of J z . There is excellent agreement with the finite-size spectrum of three-channel Kondo effect obtained from 
conformal field theory, whose energy levels (in units of the fundamental level spacing) are indicated by arrows. For J z > J z 
(right two panels), the finite-size spectrum remains Fermi-liquid-like down to minuscule energies, without crossing over to the 
non- Fermi- liquid fixed point of the three-channel Kondo effect. The persisting drift of NRG levels for J z > J z appears to be 
due to some weak numerical instability. A similar drift of levels occurs for —J Z ~3>J± > (ferromagnetic coupling). 
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FIG. 4: NRG phase boundary between the antiferromagnetic 
and the new ferromagnetic- like domain, for s — 1/2, k — 3 
and A = 5. Points on the solid line all converged to the non- 
Fermi-liquid fixed point of the three-channel Kondo effect. 
Points on the dashed line showed no such tendency up to 
TV = 120 NRG iterations, corresponding to E ~ 1CT 42 D. We 
estimate the phase boundary between the two domains to lie 
in between the solid and dashed lines. 



Eq. i|37|) . the critical coupling J* predicted by bosoniza- 
tion is equal to pJ* = 4y/3/n ps 2.2. Based on the NRG 
results for A = 3 and A = 5, we estimate the location 
of the critical coupling for a symmetric box density of 
states (i.e., for A — > 1+) to be around pJ* ~ 1.2 — 1.3. 
Thus, the bosonization and NRG results are within a 
factor of two from one another. Considering that J* lies 
well beyond the strict range of validity of bosonization, 
we find this degree of agreement to be quite remarkable. 
As for the shape of the phase boundary, it was predicted 
in Sec. IIIII to have the linear form J z = J* + Ck\ J±\, 
where < Cfc < 1 is given by Eq. As seen in Fig.0] 

the NRG phase boundary for A = 5 is well described by 
a linear curve, at least up to pj± sa 0.08. The corre- 
sponding slope Cnrg ~ 0-67 is indeed less than one, but 
is nearly three- fold larger than the bosonization result, 
Cfc = 3 = 0.25. This discrepancy can be largely accounted 
for by plugging the NRG value for pJ* into the right- 
hand side of Eq. JSSJl, which yields C w 0.5. 



We do note, however, a persisting drift of the NRG 
levels for J z > J°. A similar drift of levels occurred 
for — J z ;§> |Jj_| > (ferromagnetic coupling), and ap- 
pears to be driven by some weak numerical instability. 
While we cannot entirely rule out an eventual crossover 
to the non-Fermi-liquid fixed point of the three-channel 
Kondo effect at some lower temperature, we find this 
scenario highly unlikely given the extremely low energy 
scales reached and the rapid change of behavior as J z 
is swept through J%. We therefore identify J z with the 
phase boundary between the antiferromagnetic and the 
new ferromagnetic-like domain predicted by bosoniza- 
tion. 

Guided by this interpretation, we turned to explore the 
dependence of J z on J±. The resulting phase diagram is 
plotted in Fig. 0] for A = 5. While points on the solid 
line all converged to the non-Fermi-liquid fixed point of 
the three-channel Kondo effect, points on the dashed line 
showed no such tendency up to N = 120 NRG iterations 
(corresponding to E ~ 10 -42 -D). We therefore estimate 
the phase boundary between the antiferromagnetic and 
the ferromagnetic-like domains to lie in between the solid 
and dashed lines. In contrast to the solid line, which defi- 
nitely lies on the antiferromagnetic side of the transition, 
one cannot guarantee that all points on the dashed line 
lie on the ferromagnetic-like side. We have confirmed, 
however, that the phase diagram of Fig. 0] is practically 
unchanged upon going to N = 180 iterations. We further 
stress that the exact position of the phase boundary is 
in general A dependent, although its dependence on A 
appears to be weak. For A = 3, for example, the position 
of the phase boundary varied by no more than just a few 
percent. We expect a similar proximity of the A — > 1 + 
phase boundary. 

The above results are clearly in good qualitative agree- 
ment with the bosonization treatment of Sec. IIIII We now 
turn to a more quantitative comparison. Setting k = 3 in 



VII. DISCUSSION AND SUMMARY 

The multichannel Kondo Hamiltonian is an important 
paradigm in correlated electron systems, with possible 
applications to varied systems. Depending on the size 
of the impurity spin, s, and the number of independent 
conduction-electron channels, k, it can display either 
local Fermi-liquid, singular Fermi-liquid, or non-Fermi- 
liquid behavior. Although the isotropic model is by now 
well understood, we have shown in this paper that an 
XXZ spin-exchange anisotropy has a far more elaborate 
effect on its low-energy physics than previously appre- 
ciated. Below we briefly summarize our main findings 
and discuss their implications. A detailed account of our 
results is presented in Fig. ^ 

We begin with a spin-i impurity and with k > 2 
conduction-electron channels. From conformal field the- 
ory it is known that the non-Fermi-liquid fixed point of 
the corresponding Kondo model is stable against a small 
spin-exchange anisotropy. 15 However, it was found to be 
unstable against a sufficiently large J z ^> \ J±\ > 0. In 
the latter regime, the system flows to a line of stable 
ferromagnetic-like fixed points with a residual isospin- 
\ local moment. The phase diagram of the model thus 
consists of three distinct domains: the conventional ferro- 
magnetic and antiferromagnetic (i.e., non- Fermi-liquid) 
domains, plus a second ferromagnetic-like domain lo- 
cated deep in the antiferromagnetic regime. The new 
ferromagnetic-like domain extends above a critical lon- 
gitudinal coupling J*, whose magnitude depends on k. 
While J* exceeds the bandwidth for intermediate values 
of k, it is pushed to weak coupling for k 3> 1. Each 
of the two ferromagnetic-type domains is separated from 
the antiferromagnetic one by a Kosterliz-Thouless line, 
as depicted in Fig. ^a). While we cannot entirely rule 
out the possibility of yet another domain for sufficiently 
large | Jj_ | 3> J z , there are no indications at this point in 
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favor of such a scenario. 

Proceeding to 1/2 < s < k/2 — 1/2, spin-exchange 
anisotropy is knowni^ to be a relevant perturbation at 
the overscreened non-Fermi-liquid fixed point for fc > 4 
(for k = 4 it is a marginal perturbation). The basin 
of attraction of the overscreened fixed point is therefore 
confined to the line J z = \ J±\ > 0. The nature of the 
low-temperature fixed points for J z ^ |Jj_| was never 
explored. As shown in Sees. iHlandlVl the system flows 
to a line of stable ferromagnetic-like fixed points with 
a residual isospin-^ local moment both for sufficiently 
small J z > |Jj_| (i.e., in the limit J Z ,J± — > for any 
given ratio J Z /\J±\ > 1) and for a sufficiently large J z ^> 
\J±\ > 0. This suggests a single generic behavior in the 
entire domain J z > \J±\. Further support in favor of this 
interpretation is provided in Appendix lAl for the special 
case of a spin-one impurity. As for the domain \J Z \ < 
\J±_\, here information is confined to the weak-coupling 
regime, \pJ±_ \ <C 1/ks. Depending on the parity of 2s, the 
system flows either to a conventional Fermi liquid with no 
residual degeneracy (integer s), or to a fc-channel Kondo 
effect with an effective spin-^ local moment (half- integer 
s). It remains to be seen to what extent is this behavior 
generic to \J Z \ < \J±\. 

For isotropic antiferromagnetic exchange and k > 2, 
the two overscreened spins s = 1/2 and s = k/2 — 1/2 
share the same low-energy physics. Indeed, the corre- 
sponding Kondo models are related for J z = J± > 
through a weak-to-strong-coupling duality. It is not sur- 
prising, then, that the overscreened fixed point shows 
the same stability for both spins against a weak spin- 
exchange anisotropypi^ A similar duality appears to hold 
also for an XXZ anisotropy, at least in the range J z > 
\J±\. Indeed, for s = 1/2 the overscreened fixed point is 
stable at weak coupling (assuming J z > — | «7j_ | ) , giving 
way to a line of stable ferromagnetic-like fixed points for 
sufficiently large J z > | JjJ > 0. For s = k/2 - 1/2 the 
roles are reversed. The overscreened fixed point is stable 
against a sufficiently large J z S> | J±\ > 0, but is unstable 
for sufficiently small J z > \J±\. In the latter regime, the 
system flows to the same line of stable ferromagnetic-like 
fixed points that is approached for s = 1/2 and a large 
J z . Note that, for s = k/2 - 1/2 and small |«7j_ j > \ J Z \, 
the stability of the overscreened fixed point depends on 
the parity of k. The overscreened fixed point is stable for 
half-integer s (even fc), but is unstable for integer s (odd 
fc). 

All cases discussed above pertain to an overscreened 
impurity. We now turn to an underscreened spin, i.e., 
s > k/2. As shown in Sees. ITT1 and Ivl an XXZ anisotropy 
is a relevant perturbation both near the free-impurity 
fixed point and for a sufficiently large J z ^> |Jj_|. The 
sole exception to the rule is the case s = k/2 + 1/2, 
where 5 J = J z — |J]_| > is a marginal perturbation in 
each of these limits. Consider first the range J z > \J±\- 
Here the flow in both extremes is to the same line of 
ferromagnetic-like fixed points with a residual isospin-^ 
local moment. As before, this suggests a single generic 



behavior throughout the domain J z > |J_i_|. Similar 
to the case of isotropic antiferromagnetic exchange, the 
low-energy physics is comprised of quasiparticle excita- 
tions plus a residual local moment. However, the resid- 
ual degeneracy for s > k/2 + 1/2 is smaller than that 
of the isotropic underscreened fixed point (two versus 
2s — fc + 1 > 2), distinguishing the ferromagnetic-like 
line of fixed points from the isotropic underscreened one. 

Of the different possible cases, the most intriguing per- 
haps is that of an underscreened impurity with | J z \ < 
| J±\ and half-integer s. As we have shown in Sec. [HI the 
resulting low-energy physics is that of a fc-channel, spin- 
i Kondo effect, at least in the limit of sufficiently weak 
coupling. Thus, an XXZ anisotropy drives the system 
from underscreened to overscreened behavior. A more 
complete characterization of the transition between these 
two distinctly different behaviors is clearly needed. 

Going back to s = 1/2, we conclude with a few further 
comments on the mapping of Eq. We first reiterate 
that J* lies well beyond the strict range of validity of 
bosonization for intermediate values of fc. Nevertheless, 
this approach (and its Anderson- Yuval equivalent^) well 
describes the new ferromagnetic-like phase both qualita- 
tively and quantitatively. For fc = 3, which features the 
largest J* and is thus the most prone to error, Eq. Ij37(l 
is only a factor of two larger than the NRG estimate for 
a symmetric box density of states. We find this degree 
of quantitative agreement to be quite remarkable. 

Although the mapping of Eq. was derived in Sec. lfffl 
for a channel-isotropic model, it can easily be extended 
to channel anisotropy both in the spin-flip coupling, 
Jj_ — > Jj_ n , and in the longitudinal coupling, J z — > J zn . 
The individual transverse couplings remain unchanged in 
the course of the mapping, whether isotropic or not. As 
for the individual longitudinal couplings, these transform 
according to 

arct an (^^j - J + arctan {^£^j - 25 z , (58) 
where 

5 * = I £ arctan f^ 1 ) (59) 
n=l ^ ' 

is the average phase shift for the fc different channels. Ac- 
cordingly, the mapping of Eqs. I|58|) and l|59|l is restricted 
to values of J zn where the modulus of the right-hand side 
of Eq. H58|) does not exceed tt/2 for any of the fc chan- 
nels. Observe that channel anisotropy is preserved by 
Eqs. (|58|l and l|59|) . which adequately reduce to Eq. (JTJ 
in the limit of isotropic couplings. 

Finally, we remark on the possibility of generalizing the 
mapping of Eq. to arbitrary spin s. Two modifica- 
tions appear when the same sequence of steps is applied 
to an impurity spin larger than one-half: (i) For s > 1, 
the phase shift S z in the absence of J± is a nonlinear func- 
tion of S z . Hence, the bosonized form of the J z interac- 
tion term is no longer linear in S z as for s = 1/2. (ii) The 
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unitary transformation U produces an additional Hamil- 
tonian term of the form AS 1 ^ , similar to the one generated 
in perturbative RG [see Eq. |T2]l]. For s = 1/2, this term 
amounts to a uniform shift of the entire spectrum, which 
can be safely ignored. This, however, is no longer the 
case for s > 1/2, where different Kramers doublets are 
split. As a result of the former modification, the mapped 
Hamiltonian no longer assumes the form of a simple spin- 
exchange Hamiltonian for s > 1. The case s = 1 is an 
exception in this regard. The mapped Hamiltonian does 
acquire an additional AS*^ term, but otherwise retains 
the form of a conventional spin-exchange interaction. A 
detailed discussion of this particular case is presented in 
Appendix lAl 
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APPENDIX A: EXACT MAPPING FOR s = 1 

In this appendix, we extend the mapping of Eq. to a 
spin-one impurity. As explained in the main text, s = 1 is 
the only other spin size for which a simple spin-exchange 
interaction is restored at the conclusion of the mapping. 
However, an additional local field proportional to Si is 
generated. The basic steps of the derivation are nearly 
identical to those carried out in Sec. lIIII for s = 1/2. Only 
a few minor modifications appear, as specified below. 

Our starting point is the Hamiltonian of Eq. (|19fl . 
where S now represents a spin-one operator. Bosoniz- 
ing the fermion fields according to Eq. I|2(JI) . the bosonic 
Hamiltonian assumes the form of Eq. I|23|) with one 
sole variation: The longitudinal spin-exchange term now 
reads 



<*i2iE[ V$ «T(0)-V$ ni (0)]S z , (Al) 
71 " n=l 



replaced with 



where 



5i = arctan 



/ npJ z 
V 2 



Converting to the boson fields of Eqs. (|25|l and 



iJ f* a (0)S, 



(A2) 

and 
the 



applying the transformation U = exp 

transformed Hamiltonian Ti' = UTiW retains the same 
overall form as Eq. (|3 1 1> . but with two important mod- 
ifications: (i) The coefficient of the V$ s (0)5' z term is 



r 27T 

oi r- 



2k 



k J 2n 2 p 

(ii) A new Hamiltonian term 7Ya = AS^ with 

Si 



8D 



(A3) 



(A4) 



is added to Ti.'^ Proceeding with the transformation 
$ s (a;) — > —<S> s (x) and converting back to a fermionic 
representation, the Hamiltonian TC regains the form of 
Eq. i|19|) . but with certain renormalized parameters: 



Jz 



° Z1 



and 



A = -> 8D 



l_Si 

k 7T 



Here J' z is determined from the equation 



arctan 



irpj' z 



2vr 

T 



arctan 



irpj z 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



which comes in place of Eq. Q for a spin-i impurity. 

As for s = 1/2, the mapping of Eqs. I|A5|I -I)A9 |I is 
restricted to values of J z where the right-hand side of 
Eq. (|A9|) does not exceed 7r/2. This constrains the map- 
ping to k > 3 (overscreened impurity), and to J z > J min 
with 



Jmin ^ 







> 



k = 3 
k = 4 



(A10) 



tan(f-^)<0 k>A 



Specifically, for k > 4 the coupling regimes J m j n < J' z < 
and J z > J* > with 



2 ^ 
J, = — tan — 
irp \k 



(All) 



are mapped onto one another. This should be compared 
with Eqs. JSBJl and Egl for s = 1/2. 

Since Eqs. I)A5|1 - I|A9J| define yet another multichannel 
Kondo Hamiltonian with both spin-exchange anisotropy 
and a potentially competing AS"^ term, we have no con- 
clusive way to deduce its fixed-point structure through- 
out the J z -J± plane. Nevertheless, the behavior in one 
particular region is clear. For values of J z where both 
J' z < — \J±\ and A < 0, the spin-exchange interaction 
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FIG. 5: (Color online) Lower bound on the basin of attraction 
of the line of stable ferromagnetic-like fixed points, for s = 1 
and different values of k S> 1. In the region below and to the 
right of each solid line, both A < and J' z < — \J±\. Hence 
the system flows to the line of stable ferromagnetic-like fixed 
points with a residual degeneracy of two. The fixed-point 
structure in the remaining portion of the domain J z > \J±\, 
i.e., in between the solid lines and the overscreened line J z = 
|Jj_| (marked by a dashed line), cannot be deduced based on 
Eqs. 1A5H - 1A9H alone. 



and the AS Z term conspire to favor a ferromagnetic-like 
state with two-fold residual degeneracy. Here the resid- 



ual degeneracy originates from the S% = 1 Kramers dou- 
blet favored by A < 0. For such values of J z , one can 
safely conclude that the system flows to the line of stable 
ferromagnetic-like fixed points identified previously from 
the strong-coupling expansion of Sec.[V] These consider- 
ations provide us with the following lower bound on the 
basin of attraction of the ferromagnetic-like line of fixed 
points: 



J z > — tan 
irp 



\J±\ < Wn 



2tt {kp\J±\ 
— — h arctan — - — 
k V 2 



(A12) 



(A13) 



Figure |S] depicts that portion of the J z -J± plane de- 
fined by Eqs. (|A12|I and (|A13|) . for several values of k ^> 
1. With increasing k, a growing fraction of the domain 
Jz > \J±\ is covered by this region, which stretches for 
k 1 from pj z 4/fc <C 1 and upward in J z . This result 
partially bridges between the strong- and weak-coupling 
limits (pJ z 3> 1, \pJ±\ and < \pJ±\ < pJ z <C 1/fe, re- 
spectively), where flow to a ferromagnetic- like state has 
been established in Sees. ITU and W\ using vastly different 
techniques. As for the remaining fraction of the domain 
Jz > \J±\i its fixed-point structure cannot be immedi- 
ately deduced based on Eqs. HA5 jl -l|A9 [l alone. 
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